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ABSTRACT 

We have studied the kinematics traced by the water masers located at the centre of 
the planetary nebula (PN) K3-35, using data from previous Very Large Array (VLA) 
observations. An analysis of the spatial distribution and line-of-sight velocities of the 
maser spots allows us to identify typical patterns of a rotating and expanding ring in 
the position-velocity diagrams, according to our kinematical model. We find that the 
distribution of the masers is compatible with tracing a circular ring with a ~0'/021 
(~100 AU) radius, observed with an inclination angle with respect to the line of sight 
of 55°. We derive expansion and rotation velocities of 1.4 and 3.1 km s _1 , respectively. 
The orientation of the ring projected on the plane of the sky, at PA~158°, is almost 
orthogonal to the direction of the innermost region of the jet observed in K3-35, 
suggesting the presence of a disc or torus that may be related to the collimation of 
the outflow. 
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1 INTRODUCTION 

High angular resolution observations of molecular gas 
have revealed the presence of dense equatorial discs and 
tori towards several late asymptotic giant branch (AGB) 
stars and young planetary ne bul ae (P Ne) : for instance, 
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The interaction between the post-AGB wind 
equatorial structures has been proposed 
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and such 
as one of the 
possible physical mechanisms in shaping the bipolar and 
multipolar morphologies seen i n PN e and pro t o-PNe 
(PPNe) dBalick. Preston fc Ickei 1 19871 : iMellemal 1 19951 : 



ISahai fc Traugeij Il998l : llmai et alj bood : 
it is believed that a more effective col- 
limation mechanism(s) should be present in t he innermost 
region of the object , such as an accretion disc (|Morrislll987l : 
ISoker fc Liviol 1 1994 ISoker fc Rappaportl HoOO) or a stellar 
magnetic field produc e d by a rotating star l|Pa scoli 1985; 
IChevalier fc Luol Il994l : iGarcfa-Segural Il997l ). Furthermore, 
recent observations reveal compelling evidence for magnetic 
fields in post-AGB stars a nd PNe, associated with equato- 
rial discs and/or jets (e.g. Sabin, Ziilstra fc Greaves! 120071 : 



ISoker fc Rappaportll2000l : llckeH200a T 

The origin of molecular discs and tori around late AGB 
stars is not completely clear, b ut a possible explanation is 
the presence of a binary system (|M orri 1 19871: iLivio fc Sokeij 
1988: lTaam fc Bodenheimenl 19891 ; ISokerll2006l ). In this case, 
when one of the stars enters the AGB phase, some of the 
ejected material can be retained, generating an extended 
torus. However, in the case of the high velocity jets observed 
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IVlemmings. Diamond fc Imaill2006l ) 

The process for the formation of ac cretion discs in PPNe 
has b een investigated in detail by iReves-Ruiz fc Loped 
(1999). They find that a disc forms when a close binary 
system (with a substellar companion) undergoes common 
envelope evolution. For the case in which a low-mass sec- 
ondary is disrupted during a dynamically unstable mass 
transfer process, an accretion disc, with a radius of ~10 AU 
and a mass of ~2xl0~ 3 M^, fo r ms w ithin ~100 yr. Re- 
cently, [R^khorst^eiiemaijS (|2005h from their 3D sim- 
ulations based on a two- wind model with a warped disc, sug- 
gested that, in order to explain the observed multipolar and 
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point-symmetric shape of PNe, the required discs are quite 
small (~10— 100 AU). Furthermore, these disc-like structures 
should be dense (10 7 - 10 8 cm ) and in Keplerian rotation. 

Interferometric CO observations show larger toroidal 
molecular structures with sizes in the range of ~1000- 
6000 AU in PPNe and young PNe. These struc- 
tures seem to be systematically in expansion with a 
mean velocity of ~7 km s" 1 , such as in M 1-92 



jBuiarrabal. Alcolea fc Neril 119981). M 2-9 (|Zweigle et all 
Il997l), M 2-56 (ICastro-Carrizo et all 120021 ') . or KjPn 8 
( Forveille et al.lll998h . These velocities are comp arable to 
or be low those found in expanding AGB envelopes (|Hugginsl 
120071 ). Importantly, rotation has been observed in the Red 
Rectangle as well as slower expansion (^0.8 km s" 1 ), 
superi mposed on rotation, according to iBuiarrabal et all 
(2005). Note that the sizes of the tori are about one or 
two orders of magnitude more than the sizes of the disc- 
like structures proposed bylReves-Ruiz fc Lopez! l| 19991 ) and 
iRiikhorst. Mellema fc Ickel (|2005l )~ 

Maps of water maser emission allow the identification 
of a d isc- like structure with a radius of Q ','12 in IRAS 17347- 
3139 l|de Gregorio-Monsalvo et al.ll2004l ), whic h corresponds 
to ~ 1 00-750 AU at the source distance (|G6mez et al.l 
l2005al ). In this case, the gas kinematics suggests the pres- 
ence of both rotation and expansion in the disc traced by the 
water masers. In this context, it is very important to study 
in detail the kinematics of much smaller disc-like structures 
that might be related to the collimation of the observed 
bipolar outflows in some PNe. 

K3-35 is a young PN that sho ws a bipolar outflow in 
optical images (M iranda et al.ll20001 ). At radio wavelengths, 
K3-35 exhi bits a bright core an d two bipolar lobes with 
an S-shape l|Miranda et al.l [20011 ). The distance to t his ob- 
ject has been estimated to be ~5 kpc ( |Zhangiri995l ). using 
an statistical method. However, we note the large uncer- 
tainty of this type of estimate, since the application of dif- 
ferent statistic al methods cou ld give distances varying by 
factors of ~3 l|Phillipsl liobl ). The characteristic S-shape 
morphology of the radio lobes can be successfully repro- 
duced by a precessing jet evolv ing in a dense circumstellar 
medium ( Vela zquez et al.ll2007l ). Water maser emission has 
been found in three regions: two regions located at the tips 
of the bipolar radio jet about ~l" from the centre (regions 
N and S), and another region towards the core of the nebula 
within ^0 / .'02 (region C), suggesting the presence of a torus 
l|Miranda et al.l [200lh . In addition, OH maser emission has 
been detected towards the centre of K3-35 (within ~0'.'04), 
showing circular polarisation that suggests the presence o f 
a magnetic field (|Miranda et alJfeOQll ; iGomez et aill2005bl ). 

We decided to study the spatio -kinematical di s tribu - 
tion of the water masers reported by iMiranda et all (|200ll ) 
towards the centre of K3-35 in order to identify possible 
expansion and/or rotation motions. The paper is organised 
as follows. In Section 2, we present a simple kinematical 
model of a ring including both expansion and rotation, and 
we calculate the pattern delineated in the position-velocity 
diagrams. We also describe the least-squares fit procedure 
that we used. In Section 3, we present the current observa- 
tional data of H2O masers in the PN K3-35. We then apply 
our model to the H2O masers located towards the centre 
of the PN K3-35, making a comparison between the results 



and the observations. Finally, in Section 4, we discuss the 
implications of our results. 



2 ROTATING AND EXPANDING RING 
MODEL 

2.1 Model 

We assume a narrow, uniform, rotating, and expanding ring 
of radius R, arbitrarily oriented with respect to the line of 
sight. Its projection on the plane of the sky is an ellipse 
with semimajor and semiminor axes a and 6, respectively. 
We define the two frames of reference shown in Fig. 1. Both 
coincide with the plane of the sky, one of them has its origin 
at the centre of the ellipse and is oriented such that the x'- 
axis is along the major axis of the projected ellipse, and the 
other one has the axes parallel to the RA and Dec axes. The 
semimajor and semiminor axes are related to the ring radius 
by a = R and b — Rcosi, where i is the inclination angle 
between the line of sight and the normal to the ring plane 
as shown in Fig. 2. 

The equation of the ellipse is given by 
12 12 

^ + V = 1 > w 

and the transformation equations between the coordinate 
systems are 



x = (x — xo) cos 6 + (y — yo) sin 8, 
y' = —(x — x ) sin 6* + (y - y ) cos* 



(2) 
(•3) 



where (xo,yo) is the position of the centre of the ellipse, 
and 6 is the angle between the :r-axes of the two frames of 
reference, and is defined as positive clockwise. The angle 9 
is related to the position angle (PA) of the major axis of the 
ellipse by PA = 90° - 6. 

Let v„ , v Iot , and w cxp be the local standard of rest (LSR) 
systemic velocity, the rotation velocity, and the expansion 
velocity of the ring, respectively. Then the observed LSR 
velocity of a point in the ring can be expressed as 



SR 



v s H «rot sini + 



y 



, tan 1. 



(4) 



Hence the observed Vlsr will be a linear function of either 
x or y' , if only one type of motion ( rotation or expansion , 
respectively) is present in the ring l|Uscanga et al.l [20051 ). 
Using equation (1), equation (4) can be written either in 
terms of the x' or y' coordinate as 



[Vlsr. - v* - (x'/a)v I0 t sini]' 
(vexp sini) 2 

[Vlsr. — v s — (y'/a)v exp tani] 



a 2 



+ 



V ' 



(5) 



(6) 



(u rot sini) 2 (ocosi) 2 

Therefore, equations (5) and (6) indicate that the observed 
Vlsr has a quadratic form (ellipse) expressed in terms of x' 
or y' , when both motions are present. 

In this model, we do not solve the radiative transfer 
through the ring, but we assume that the emission at a given 
Vlsr comes from the point of the ring having this line-of- 
sight velocity component. We then use this information to 
construct position-velocity diagrams. 

The orientation of the major axis of the ellipse in the 
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Figure 2. Position- velocity diagrams for a rotating and expanding ring model. The top panels correspond to a positive inclination angle 
and the bottom panels correspond to a negative inclination angle. The sense of rotation (clockwise or counterclockwise) as seen from the 
observer is indicated. The x' and y' axes are those defined in Fig. 1. 
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Figure 1. Reference systems. Both the x'y'- and the xy- 
coordinate systems are in the plane of the sky. The x'y'-system 
has its origin at the centre of the ellipse (xo, yo), and 8 is the 
angle between the x-axis and the x'-axis. The x and y axes are 
parallel to the RA and Dec axes. 

position-velocity (x'-Vlsr or j/'-Vlsr) diagrams changes de- 
pending on the value of the inclination angle and on whether 
the sense of rotation is clockwise or counterclockwise as seen 
from the observer's point of view, as shown in Fig. 2. Accord- 
ingly, we are able to distinguish between a positive or neg- 
ative value of the inclination angle and the sense of the ro- 
tation by doing a comparison between the position-velocity 
diagram delineated by the maser emission and the position- 



velocity diagrams expected for a rotating and expanding 
ring. It is important to note that similar position-velocity 
diagrams can be obtained considering contraction instead 
of expansion in the ring, but changing the sign of the incli- 
nation angle and the sense of rotation. This ambiguity can 
be solved by constraining the value of the inclination an- 
gle (positive or negative) with additional information (see 
Section 3.2). 

2.2 Least-squares fit 

We carried out a least-squares fit of an ellipse to the observed 
emission, to estimate its spatial distribution on the sky. In 
order to do this, we considered the curve on the a;y-plane 
given by equations (l)-(3) for a given set of values of the 
parameters (xo,yo) (the centre of the ellipse on the plane 
of the sky), a and b (the semimajor and semiminor axes, 
respectively) and 9 (the angle between the x-axes of the two 
frames of reference). We then compute the minimum dis- 
tances dj between the position (xj,yj) of each of the masers 
and the ellipse (these distances are measured along straight 
lines that pass through the maser and intersect the ellipse at 
right angles to the curve). With these distances, we define 
the x 2 as 

N / \ 2 

• 

where N is the number of data points (i.e. the factor N — 5 
is the number of degrees of freedom), and a-j is the error 
associated with the position of the maser spots. We then 
find the set of parameters (xo,yo,a,b,0) (see above) which 
give the minimum of \ 2 ■ 

As a first approximation, we began by setting the values 
of aj equal to the observational uncertainties Aj for the mea- 
sured positions of the masers, and then finding the ellipse 
that gives the minimum value of \ 2 ( see equation [7}. The 
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actual structure of the maser-emitting region could be a ring 
of finite width, for which a broad ellipse is a rough approx- 
imation of its projection on the plane of the sky. Therefore, 
we do not expect maser spots to exactly trace an ellipse. We 
can characterise the width of the ring, assuming that the fit- 
ted ellipse traces the mean projected angular distance of the 
maser emission to the central star, and the actual emission 
will be distributed around this ellipse, with a dispersion A e . 
We treat A e as a source of error for the ellipse fit, additional 
to the measured error, so that aj 2 = Aj 2 + A e 2 . We then 
try different values of the width parameter (A e ), until we 
obtain a fit with a minimum x 2 (A e ) = 1. We consider 2A e 
as the characteristic width of the maser ring. 

The (xo,yo,a,b,0) parameters obtained from the min- 
imization of \ 2 gi ye the best elliptical fit to the observed 
positions of the masers. With these spatial parameters, we 
carried out a kinematical fit, using the LSR velocities of the 
maser components in order to define a^ 2 , of the form 

JV 

2 1 V- 1 

X V 



0.03 



TV - 3 ^ a v 2 

3=1 



x Vlsr, j - «s 



-Vrat SUl J 



(8) 



where TV — 3 is the number of degrees of freedom, and 
(x'j,y' a) are given by equations (2)-(3). Here a v is the un- 
certainty in the observed LSR velocity that we adopt as the 
spectral resolution of the observations. The minimization of 
X 2 V yielded the best values for the systemic (v s ), rotation 
(frot), and expansion (w cxp ) velocities. 



3 H z O MASERS IN K3-35 (REGION C) 
3.1 Observational Data 

There are only two sets of VLA water maser observations 
towards the PN K3-35. The VLA 1 999.7 epoch observa- 
tions reported bv lMiranda et all (l200ll), and the VLA 2002.3 
epoch observations reported bv lde Gregorio-Monsalvo et al.l 
(2004). In the latter observations, only a group of four maser 
spots was detected towards the central region of this source. 
No maser emission was detected at the tips of the bipolar 
lobes of the PN. 

Note that in the lde Gregorio-Monsalvo et al.l (|2004T ) pa- 
per the position of the continuum peak used to align the 
positions of the masers at the two epochs was not used with 
enough precision, resulting in a spurious shift of the maser 
spots from one epoch to another (see their Fig. 4). Using the 
position with the adequate precision, we find the positions 
of the masers at the two epochs to be consistent within the 
uncertainties, O'.'Ol (2a). 



3.2 Model Application and Results 

In our analysis, we have used the water maser data from the 
VLA 1999.7 epoch o bservations tow a rds th e central region 
of K3-35 reported bv lMiranda et ail l|200ll ). At this epoch, 
the number of maser spots detected was larger than dur- 
ing the other epoch, allowing us a better identification of 
possible expansion and/or rotation motions at the centre of 
K3-35. The velocity resolution of the VLA observations was 
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Figure 3. Positions of the K3-35 water maser spots in offsets 
relative to the positi on of the 1.3 cm continuum emission peak 
iMiranda et alj|200ll) . Each spot is labelled with its correspond- 
ing LSR velocity (in km s — 1 ). The dashed ellipse corresponds 
to the least-squares fit to the maser spots positions whose pa- 
rameters are indicated in Table 2. The open circle indicates the 
nominal position of the 1.3 cm continuum emission peak (see Ta- 
ble 1), its size is equal to the uncertainty of this position. The 
diamond indicates the position of the centre of the ellipse that 
was obtained from the fit. The straight line shows the direction 
of the bipolar outflow traced by the innermost region of the jet. 
The arrows show the sense of rotation of the proposed ring (see 
Section 2.1). The broad cyan ellipse has a width 2A e , with A e 
fulfilling x 2 (A e ) = 1 (see Section 2.2). The positional error bars 
indicate the uncertainties in the relative positions between maser 
spots given in Table 1. 



1.2 km s _1 and the accuracy in the relative positions of the 
water maser spots was of the order of milliarcseconds. The 
positions of the observed water maser spots towards the core 
(region C) are listed in Table 1. We adopt as the origin of the 
xj/-coordinate system the position of the 1.3 cm continuum 
emission peak. 

Based on a least-squares fit to the positions of the maser 
spots (see Section 2.2) and using the observational uncer- 
tainties given in Table 1, we have found that the H2O masers 
located towards the core of K3-35 can be fitted by a circular 
ring of R ~ 0'.'021 radius (~100 AU at the estimated dis- 
tance of ~5 kpc) with an angular width of 2A e = 0'.'003, 
observed at an inclination angle of \i\ ~ 55° (see Table 2 
and Fig. 3). 

Spectroscopic observations show that the northeastern 
lobe of the ou tflow is blueshifted, a nd the southwestern one 
is redshifted (|Miranda et alj feoooh. If we assume that the 
ring traced by the water masers is perpendicular to the 
bipolar lobes then the inclination angle should be positive 
(i ~ +55°). This means that the western half of the ring is 
closer to the observer. 

The calculated position of the centre of the ellipse rela- 
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lable ±. 


H2O Masers in 


K3-35 (Region C) 






VlSR 


Flux Density 


Position 


Position Uncertainty^ 


(km s _1 ) 


(mjy) 


a(J2000) 


<5(J2000) 


(arcsec) 


24.6 


23 


19 27 44.0243 


21 30 03.438 


0.010 


24.0 


61 


19 27 44.0242 


21 30 03.428 


0.004 


23.3 


218 


19 27 44.02246 


21 30 03.4460 


0.0012 


22.6 


1010 


19 27 44.022254 


21 30 03.45146 


0.00024 


22.0 


1572 


19 27 44.022364 


21 30 03.45335 


0.00014 


21.3 


945 


19 27 44.02247 


21 30 03.4564 


0.0003 


20.7 


201 


19 27 44.02257 


21 30 03.4625 


0.0011 



a Units of right ascension are hours, minutes, and seconds, and units of declinat ion are 
degree s, arcminutes, and arcseconds. Data from lMiranda et al.l j200lh and lGomez et all 
l|2003h . 

" Relative position uncertainties (2a) between maser spots. The position of the 1.3 cm 
continuum emission peak is a(J2000) = 19 h 27 m 44?0233, <5(J2000) = 21°30'03'/441. 
The relative position uncertainty between the continuum and the H2O masers is 0"002. 
The accuracy of the absolute positions is 0'.'05. 



tive to the position of the 1.3 cm continuum emission peak 
(xo,yo) is given in Table 2. Both positions are in agreement 
within the uncertainties (see the overlap of these positions 
in Fig. 3). 

The kinematic trend is shown in Fig. 4. Since we have 
determined that the inclination angle of the ring is positive, 
the ambiguity between expansion and contraction can be 
solved when we compare the position-velocity diagrams de- 
lineated by the water masers (see the top panel of Fig. 4) and 
the position- velocity diagrams of the model (see Fig. 2). We 
have found that the ring traced by the masers rotates clock- 
wise as seen from the observer at a velocity v TO t — 3.1 km s _1 



and expands at a velocity v e 



1.4 km s" 1 (see Table 2). 



The rotation and expansion velocities are estimated from 
a purely kinematical fit to the LSR velocities of the maser 
spots (see equation 8). The kinematical fit yields a x 2 v value 
of 1.94, which means that the fit is good, assuming a conser- 
vative confidence level of 90%. Although not all maser spots 
may be completely independent, given the limited angular 
and spectral resolution, the observed velocity trend and the 
good kinematical fit suggest that these motions are real and 
systematic. We note that we have considered a single value 
of v rot instead a Keplerian rotation law. This is reasonable, 
since the velocity gradient over the ring width would be only 
~7%. Tracing more subtle velocity variations would require 
more data points than the ones available. 

The expansion velocity we found for the ring is close 
to that of thermal motions or subsonic turbulence. However 
we do not expect that thermal or turbulent motions could 
produce the systematic motions described in the position- 
velocity diagrams. The water masers seem to be tracing 
a spatio-kinematical structure with organised motions (at 
macroscopic scales). 



4 DISCUSSION 

From our model, we conclude that a ring is a likely explana- 
tion for the distribution and kinematics shown by the water 
masers located towards the centre of the PN K3-35. This 
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Figure 4. Position-velocity diagrams. Top: The ordinate axis 
corresponds to the observed LSR velocity and the abscissa axis 
corresponds to the coordinate x' or y' . Bottom: Same as Fig. 4 
(top), but the abscissa axis corresponds to right ascension offset or 
declination offset relative to the position of the 1.3 cm continuum 
emission peak. The points correspond to the observed maser spots 
towards the centre of K3-35 and the dashed ellipses correspond to 
the kinematical model using the parameters listed in Table 2. The 
error bars in position are those shown in Fig. 3. The uncertainty 
in velocity is ~0.6 km a -1 . 



ring may be arising from the innermost region of a disc or 
torus probably formed at the end of the AGB phase. Since 
masers trace regions with very stringent physical conditions, 
it is not possible to know whether the ring is tracing part 
of a toroidal or a disc-like structure. The kinematics of the 
ring in K3-35 suggests the presence of both rotating and 
expanding motions, as was also proposed in the young PN 
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Table 2. Results of the rotating and expanding ring model fit to 
the H2O masers towards the centre of K3-35. 



Parameter Value" 



a 


0"021 


± 


07003 


b 


07012 


± 


07002 


x * 


07001 


± 


07001 


yo * 


07004 


± 


07004 


PA 


158° 


± 


10° 


i 


55° 


± 


7° 


Vs 


22.8 


± 


0.5 km s — 1 


^cxp 


1.4 


± 


0.9 km s" 1 


ffot 


3.1 


± 


0.8 km s" 1 



a Note: Uncertainties are 2a. 

* Coordinates of the centre of the ellipse relative to the position 
of the 1.3 cm continuum emission peak. 

IRAS 17347-3139 l|de Gregorio-Monsalvo et al.ll2004T ). The 
estimated expansion and rotation velocity values for K3-35 
are similar (a few km s _1 ) to those obtained from water 
maser observations of IRAS 17347-3139. 

The calculated expansion velocity of the ring 
(~1.4 km s _1 ) in K3-35 is comparable to, but below, 
the expansion velocities of the tori inferred from interfer- 
ometric CO observations in some PPNe and young PNe. 
For instance, the expansion velocities in M 1-92, M 2-9, 
M 2-56, a nd KjPn 8, have values in the range of ~5- 
8 kin s' 1 llBuiarrabal. Alcolea fe Nerilll99Sl; IZweigle etafl 
ll997l ; ICastro-Carrizo et al.ll2002l ; iForveille et alJll998T ). litis 
noteworthy that the estimated expansion velocity in K3-35 
obtained from water maser observations is quite similar to 
the slow expansion velocity (~0.8 km s _1 ) deduced from CO 
obser vations in the PPN Red Rectangle (jBuiarrabal et ail 
l2005h . However, the estimated sizes of the tori detected in 
CO observations in those sources are about one order of 
magnitude more than the size of structure traced by the 
water masers in K3-35, suggesting that the structures ob- 
served in CO are probably related to the outermost equa- 
torial region, while the water masers could be arising from 
the innermost one, such as a circumstellar disc, closer to the 
central star. 

In the case of K3-35, the water masers may be delineat- 
ing an ellipse on the plane of the sky with a semimajor axis 
of ~100 AU and a PA~158±10°, suggesting the presence of 
a disc or torus projected on the sky. The major axis of the 
ellipse is almost perpendicular (within the uncertainties) to 
the direction of the bipolar emission traced b y the innermost 
region of the jet, with a PA~65°observed by I Miranda et all 
l|200ll ). suggesting that the disc or torus could be physically 
related to the collimati o n of t he bipolar outflow. 

Recently, iHugginsl (120071 ) found that jets typically ap- 
pear a few hundred years after the torus formation in PPNe 
and young PNe. This time sequence provides evidence that 
jets and tori ar e phys ically related. In the case of K3-35, 
IVelazquez et al.l (|2007T l modelled its radio continuum emis- 
sion as a precessing dense jet with an age of ~40 yrs. On 
the other hand, our estimate of the dynamical age of the 
ring traced by the maser emission is about ^350 years, us- 
ing its radius (~100 AU) and assuming a uniform expansion 
velocity of ~1.4 km s _1 . These values would, in principle, 
suggest a lag in time similar to the value found by IHugginsl 



However, the dynamical age we derive is not reliable 
enough as an age estimate. Proper motion studies of the wa- 
ter maser emission (e.g. using e-MERLIN and VLBA) could 
provide better estimates. 

The kinematics of the ring in K3-35 suggests the pres- 
ence of both expansion and rotation. A rough estimate for 
the central stellar mass can be obtained by assuming that 
the total energy (kinetic plus gravitational) of the masing 
gas is close to zero. In this case, M ~ (R/2G)(v^ p + v^ ot ), 
where R ~ 100 (D/5 kpc) AU, is the radius of the ring. 
Hence, M ~ [0.7 (D/5 kpc) ± 0.3] M , where D is the 
source distance. The error in the mass only includes the er- 
rors in the fitted parameters. This estimate of the central 
mass is in agreement with the core mass requir ed for a PNe, 
according to evolutionary models (Kwok 2003). 



ACKNOWLEDGMENTS 

LU and YG acknowledge support from DGAPA-UNAM 
grant IN100407 and CONACyT grant 49947. LU is sup- 
ported by Secretan'a de Estado de Universidades e Inves- 
tigation of MEC. AR and JC are supported by CONA- 
CyT grants 46828-F and 61547. LFM acknowledges support 
from MEC AYA2005-01495 grant (co-funded with FEDER 
funds). GA, JFG, and JMT are supported by the MEC 
AYA2005-05823-C03 grant (co-funded with FEDER funds). 
GA, JFG, LFM, JMT and LU are also supported by Con- 
sejeria de Innovation, Ciencia y Empresa of Junta de An- 
daluda. We thank L. F. Rodriguez for valuable comments. 
We are thankful to our referee, Anita Richards, for her useful 
comments on the manuscript. 



REFERENCES 

Balick B., Preston H. L., Icke V., 1987, AJ, 94, 1641 
Bieging J. H., Nguyen-Quang-Rieu, 1988, ApJ, 324, 516 
Bujarrabal V., Alcolea J., Neri R., 1998, ApJ, 504, 915 
Bujarrabal V., Castro-Carrizo A., Alcolea J., Neri R., 2005, 

A&A, 441, 1031 
Bujarrabal V., Neri R., Alcolea J., Kahane O, 2003, A&A, 

409, 573 

Castro-Carrizo A., Bujarrabal V., Sanchez Contreras O, 

Alcolea J., Neri R., 2002, A&A, 386, 633 
Chevalier R. A., Luo D., 1994, ApJ, 421, 225 
de Gregorio-Monsalvo I., Gomez Y., Anglada G., Cesaroni 

R., Miranda L. F., Gomez J. F., Torrelles J. M., 2004, 

ApJ, 601, 921 
Feibelman W. A., 1985, AJ, 90, 2550 

Forveille T., Huggins P. J., Bachiller R., Cox P., 1998, ApJ, 
495, Llll 

Garcfa-Segura G., 1997, ApJ, 489, L189 

Gomez J. F., de Gregorio-Monsalvo I., Lovell J. E. J., 

Anglada G., Miranda L. F., Suarez O., Torrelles J. M., 

Gomez Y., 2005a, MNRAS, 364, 738 
Gomez Y., Miranda L. F., Anglada G., Torrelles J. M., 

2003 in Kwok S., Dopita M., Sutherland R., eds, Proc. 

IAU Symp. 209, Planetary Nebulae: Their Evolution and 

Role in the Universe. Astron. Soc. Pac, p. 263 
Gomez Y., Tafoya D, Anglada G., Franco-Hernandez R., 



7 

Torrelles J. M., Miranda L. F., 2005b, Mem. Soc. Astron. 

Italiana, 76, 472 
Huggins P. J., 2007, ApJ, 663, 342 
Icke V., 2003, A&A, 405, Lll 

Imai H., Obara K., Diamond P. J., Omodaka T., Sasao T., 

2002, Nat, 417, 829 
Kwok S., 2003 in Kwok S., Dopita M., Sutherland R., eds, 

Proc. IAU Symp. 209, Planetary Nebulae: Their Evolution 

and Role in the Universe. Astron. Soc. Pac, p. 3 
Livio M., Soker N., 1988, ApJ, 329, 764 
Mellema G., 1995, MNRAS, 277, 173 

Miranda L. F., Fernandez M., Alcala J. M., Guerrero M. 

A. , Anglada G., Gomez Y., Torrelles J. M., Aaquist O. 

B. , 2000, MNRAS, 311, 748 

Miranda L. F., Gomez Y., Anglada G., Torrelles J. M., 

2001, Nat, 414, 284 
Morris M., 1987, PASP, 99, 1115 
Pascoli G., 1985, A&A, 147, 257 
Phillips J. P., 2004, MNRAS, 353, 589 
Reyes-Ruiz M., Lopez J. A., 1999, ApJ, 524, 952 
Riera A., Garcfa-Lario P., Manchado A., Bobrowsky M., 

Estalella R., 2003, A&A, 401, 1039 
Rijkhorst E.-J., Mellema G., Icke, V., 2005, A&A, 444, 849 
Sabin L., Zijlstra A. A., Greaves J. S., 2007, MNRAS, 376, 

378 

Sahai R., Trauger J. T., 1998, AJ, 116, 1357 

Sahai R., Hines D. C, Kastner J. H., Weintraub D. A., 

Trauger J. T., Rieke M. J., Thompson R. I., Schneider 

G., 1998, ApJ, 492, L163 
Soker N., 2006, ApJ, 645, L57 
Soker N., Livio M., 1994, ApJ, 421, 219 
Soker N., Rappaport S., 2000, ApJ, 538, 241 
Taam R. E., Bodenheimer P., 1989, ApJ, 337, 849 
Uscanga L., Canto J., Curiel S., Anglada G., Torrelles J. 

M., Patel N. A., Gomez J. F., Raga A. C, 2005, ApJ, 634, 

468 

Velazquez P. F., Gomez Y., Esquivel A., Raga A. C, 2007, 

MNRAS, 382, 1965 
Vlemmings W. H. T., Diamond P. J., Imai H., 2006, Nat, 

440, 58 

Zhang C. Y., 1995, ApJS, 98, 659 

Zweigle J., Neri R., Bachiller R., Bujarrabal V., Growing 
M., 1997, A&A, 324, 624 

This paper has been typeset from a TjrjX/ L^IgX file prepared 
by the author. 



